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We present a study of the inelastic decay rate of quasiparticles in a two-dimensional electron 
gas with spin-orbit interaction. The study is done within the G°W° approximation. The spin- 
orbit interaction is taken in the most general form that includes both Rashba and Dresselhaus 
contributions linear in magnitude of the electron 2D momentum. Spin-orbit interaction effect on 
the inelastic decay rate is examined at different parameters characterizing the electron gas and the 
spin-orbit interaction strength in it. 
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I. INTRODUCTION 

Nowadays, in condensed matter physics and semicon- 
ductor microelectronics, two-dimensional (2D) electron 
system is one of the main objects of detailed study. Such 
a system is formed by, e.g., surface-state electrons or elec- 
trons in semiconductor heterostructures. Phenomenon 
that is observed in such systems and makes them of 
great interest, especially in context of spintronic appli- 
cations, is spin-orbit interaction (SOI). This interaction 
arise from the structure inversion asymmetry of poten- 
tial confining the electron system in directions perpen- 
dicular to the confinement plane (the Rashba spin-orbit 
interaction^) and the bulk inversion asymmetry that is 
present in semiconductor heterostructures based on ma- 
terials with a zinc-blende structure (the Dresselhaus spin- 
orbit interaction 2 -! 3 -) . The Dresselhaus interaction de- 
pends on semiconductor material and growth geometry, 
whereas the interaction strength of the Rashba SOI can 
be tuned via an externally applied electric field perpen- 
dicular to the confinement planed As a result, one can 
controllably manipulate the spin in devices without re- 
course to an external magnetic field 

In order to efficiently exploit the mentioned phe- 
nomenon, a theoretical study of dynamics of electrons 
and holes in the 2D spin-orbit coupled electron sys- 
tems is needed. The most discussed and studied pro- 
cesses concerning this problem are spin relaxation and 
spin dephasingji However, to our knowledge, such crucial 
quasiparticlc property as the lifetime caused by inelas- 
tic electron-electron scattering remains still insufficiently 
studied. To all appearance the first attempt to analyze 
what effect the SOI has on the quasiparticlc lifetime has 
been made in Ref. H]. In the work cited, a particular 
case of the 2D electron gas (2DEG) with the Rashba 
SOI was considered at the limit of En <C Ep, where Ep 
is the Fermi energy and En = m*a 2 /2 with a and m* 
being the interaction strength and the effective electron 



mass, respectively (unless stated otherwise, atomic units 
are used throughout, i.e., e 2 = h = m = 1.). Within 
the G°W° approximation, it has been shown that in a 
small vicinity of Ep a modification of the lifetime due 
to the SOI is insignificant and does not depend on the 
subband index of the spin-orbit split band. To go beyond 
the limits of Ref. [H, in Ref. the inelastic lifetime (de- 
cay rate) of quasiparticles in the 2DEG with the Rashba 
SOI has been studied within a wide energy region. For 
material parameters typical for In x Gai_ x As 2DEGs, it 
has been revealed that modifications induced by the SOI 
and the dependence on the subband index become no- 
ticeable, when the decay channel due to plasmon emis- 
sion appears. The first joint theoretical and experimental 
investigation of hole lifetimes in a 2D spin-orbit coupled 
electron system has been done in Ref. [1(1 In addition 
to a demonstration of the weak influence of the SOI on 
hole lifetimes by the case of the Au(lll) surface state, a 
hypothetical system, where the SOI can have a profound 
effect, has been considered. 

In this work, we generalize the results on effect of the 
SOI on the quasiparticlc lifetime. Within the G a W° 
approach with the screened interaction W° evaluated 
in the random phase approximation (RPA), we study 
the inelastic decay rate of quasiparticles in a 2DEG 
with the Rashba and Dresselhaus interactions linear in 
k — magnitude of the electron 2D momentum k. In 
our G°H /0 -calculations, material parameters suitable for 
InAs quantum wells are taken. We compare the inelas- 
tic decay rates calculated at different ratios between the 
interaction strengths of the mentioned spin-orbit inter- 
actions. We show that on the energy scale, for the taken 
material parameters, the main visible effect induced by 
the SOI is modifications of the plasmon-emission decay 
channel via the extension of the Landau damping region. 
We also consider a hypothetical small-density case, when 
in the 2D spin-orbit coupled electron system the Fermi 
level is close to the band energy at k = 0. For such 
a system, we predict strong subband-index dependence 



2 



n^-isin^t+jcos^ 




FIG. 1: (Color online) A spatial location of the vectors Uk and 
nk determining the spin-quantization axis with polar angles 
0k and $k and the rotation axis, respectively. 
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and anisotropy of the inelastic decay rate for electrons 
and appearance of a plasmon decay channel for holes. 



II. APPROXIMATIONS 

We consider a 2DEG described by the Hamiltonian 
H = Hq + Hso with Hq — k 2 /2m* and the spin-orbit 
contribution 



Hso = a (a x k y 



Oyk x ) + P (a x k x 



CTyky) 



(1) 



that includes both Rashba and Dresselhaus terms. The 
latter is written with the assumption that a quantum 
well grown in [001] direction is considered. In Eq. ([TJ), 
k x , y are the electron momenta along the [100] and [010] 
cubic axes of the crystal, respectively, a XiV are the Pauli 
matrices, m* is the effective electron mass, a and P are 
the interaction strengths for the Rashba and Dresselhaus 
spin-orbit interactions. To bring the Hamiltonian to a di- 
agonal form, we perform the rotation in spin space gen- 
erated by t/k = cxp[i(er • rik)Ok/2] dependent on the 
momentum k. The rotation is performed with the angle 
0k around the axis determined by nk- A positional re- 
lationship of the axis nk and the spin-quantization axis 
Uk is shown in Fig. [TJ We suppose that we deal with the 
in-plane spin polarization, i.e., Ok = 7r/2. In the new, 
unitary transformed, spin basis the spin-orbit contribu- 
tion has the formii 



H 'so 



UlHsoU* 

= -k [a sin(^ k - $k) + P cos((^ k + $k)] or*, (2) 

where the angle $k is related to the polar angle ip^ of 
the momentum k as 



tan$ k 



a cos (fk + P sin tp^ 
a sin ifk + P cos ip^ 



(3) 



FIG. 2: (Color online) Fermi contours and spin orientations 
in the momentum plane for different values of interaction 
strengths of the Rashba and Dresselhaus spin-orbit interac- 
tions. 



Due to the diagonal form of H' SOl the energy bands are 
simply given by^ 



2m* 



sk [a sin((^ k - $k) + P cos(^ k + $k)] (4) 



,ikr| 



and correspond to the wave functions V'ksl 1 ') 
s) with the subband index s = ±(4-, t); where I, j" arc 
the spin components in the new spin basis. This means 
that for the initial, untransformed, Hamiltonian we have 
the following eigenstates ^ks(r) = U^ip^^) . The spin 
orientation in k space reads as (see Fig. [2]) 



(V>ks|o-|V>k s ) = 




(5) 



Note that the case with a^0 and P = (pure Rashba) 
is characterized by the angle $k = <^k — tt/2, whereas in 
the situation with P ^ and a — (pure Dresselhaus) 
one has $k = 27r — ip^. In the special case of a = P, the 
angle $k = — tt/4 does not depend on tp±. 

The inelastic decay rate (inverse lifetime r s _1 (k) 
caused by inelastic electron-electron scattering) is de- 
termined by the imaginary part of the matrix ele- 
ments of the quasiparticle self-energy (S s (k, lj)) = 
(V'k s ( r i)l s ( r i) r 2;w)|'!/'k s (r2))r 1 r 2 at the energy u = E^ s 
as r s (k) = 2|Im(S s (k,£ , ks ))|. At the Hartree-Fock (HF) 
mean-field level, these elements are totally real and have 
the form 



:(k-q), (6) 
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where v c (k) = 27r/(|k|eo) is the bare Coulomb interaction 
with e being the static dielectric constant. The factors 
F K P = [! + ss ' u k • Up] /2 come from | (s'|[/t[/ k |s) | 2 and 
/ks is the Fermi factor. Such a form ([5]) is similar to 
the exchange contribution to the single-particle energies 
considered in Refs. H^ - Hil in the pure Rashba case. How- 
ever, in order to examine quasiparticle lifetimes, one has 
to go beyond the HF approximation. The simplest vari- 
ant is the G°W° approximation (for details about the 
approximation, we refer the reader to Refs. and [13). 
Within such an approximation, we arrive at the follow- 
ing expression for the imaginary part of the mentioned 
matrix elements 

Im(E s (k, W )) = - W Jj£/ qa , (7) 

s' 

x ImW (k - q,uj - E^eiE^, - u), 
when tjj < Ep, and 

Im(S s (k, W )) = J2j (02 F k S q [l-.V] (8) 

x haW°{k-q,u-E^)e(u-E^), 

when cu > Ep. In these equations, 6{x) is the step func- 
tion and the screened interaction 



^°(q, W )=« c (q)[l-P°(q, W H(q)] 1 



(9) 



is defined by the RPA irreducible polarizability (see also 
Refs. [H and [H, where the retarded part of P° was ex- 
amined) 



P°(q,a,) = E/ 



dk 



(2tt) 2 k ' k+q 

(1 - /k+qs)/ks 



lo + E ks * - -E k +q.s + ifl 

/k+qs(l - fks') 



(10) 



E 



k+qs — IT) 



III. RESULTS AND DISCUSSION 

Using the example of an InAs quantum well, we take 
the effective mass m* = 0.023 (see, e.g., Ref. l2(i) and the 
static dielectric constant e = 14.55 (see, e.g., Ref. I21I ). 
The interaction strength of the Dresselhaus SOI is cho- 
sen to be j3 = 1.6 x 10 -11 eVm to simulate a quite 
narrow quantum well^S At the ratio a/f3 = 2.4 (see, 
e.g., Ref. l22t ). we have the Rashba interaction strength 
a = 3.8 x 10~ n eVm. The electron density is put at 
112D = 2.55 x 10 11 cm -2 that corresponds to the Fermi 
energy Ep ss 26 meV. 

Fig. [3] shows our results on the inelastic decay rate 
r s (k) obtained with the material parameters listed 
above. Two main points caused by the spin-orbit split- 
ting of the band make the considered 2DEG different 
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FIG. 3: (Color online) Inelastic decay rate T s as a function 
of k at several values of the polar angle ip]< in the case of the 
ratio a//3 — 2.4. Inset: The corresponding energy bands E^ 3 
(at the left) and the same T s as a function of Ek s measured 
from the Fermi energy (at the right). Also, for reference, the 
case of the 2DEG without the SOI (a — j3 = 0) is presented. 



from that without the SOI. These arc an angle-dependent 
relative shift of V + and T_ on the momentum scale and 
some smoothing of sharp forms of the peak caused by 
opening of the plasmon decay channel. The former re- 
flects the fact that subbands of the split band reach the 
same energy at different momenta, while the latter orig- 
inates from the extension of the Landau damping region 
(the region where plasmons decay into single-particle 
excitations 2 ^) due to appearance of inter-subband transi- 
tions (for a detailed discussion of the screening properties 
of the 2DEG with the SOI we refer the reader to Refs. [U 
and fl9| ) . This extension varying with the polar angle y>k 
leads to a nonzero plasmon linewidth, when the plasmon 
spectrum enters into the SOI-induced damping region. 

In order to show what effect the SOI has on the in- 
clastic decay rate for different subbands, in the inset of 
Fig. |31 by setting up a correspondence between T s (k) 
and E\l S via the momentum k, we plot the decay rate as 
a function of energy. On first glance, it may seem that we 
have an ordinary energy dependence of the decay rate as 
in a 2DEG without the SOI: the quadratic behavior with 
the logarithmic enhancement in the vicinity of the Fermi 
energy with F s = at Ep and the jump above the Fermi 
energy, which is caused by opening the plasmon decay 
channel for excited electrons^! However, on examining 
the energy dependence of T s in detail, we can say that 
due to the finite plasmon linewidth the plasmon decay 
channel manifests itself at lower energies, when it occurs 
in a 2DEG without the SOI. The same reason leads to 
reduction in the jump. Also, we can reveal distinctions 
between T + and r_, which become noticeable, when the 
plasmon-cmission decay channel appears, and increase 
upon moving from (p^ — 7r/4 to 37r/4. An analysis of the 
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FIG. 4: (Color online) Same as in Fig. O but for the ratio 
a//3 — 1. Also, the pure Rashba and Dresselhaus cases are 
presented. 



inelastic mean-free path (IMFP) A s (k) = |V k -Ek s |/r s (k) 
as a function of energy has shown that, as a consequence 
of the distinctions between T + and r_, the IMFP of elec- 
trons can vary with the subband index. For example, at 
= 37r/4 for electrons this variation can reach, e.g., 
~ 8%. 

The obtained results can be understood by inspecting 
constant-energy contours shown in Fig. [2] with mental 
drawing of possible transitions selected by the factors 
^k S q °f Eqs. (|7|) and (J5J at a given ip^ (see also Ref.Ql). 
Actually, for each subband (s = ±) one has a set of intra- 
and inter-subband transition momenta as arguments of 
ImFT . For the chosen material parameters and for ip^ = 
7r/4, these momenta do not vary considerably with the 
subband index s. For ip-k = Stt/A differences in both 
intra- and inter-subband transitions for s = + and s = — 
become already sensible for values of ImVF , especially in 
the vicinity of plasmon peaks of the latter. 

Now, remaining m* , eq, ri2D, and Ep unchanged, we 
consider the case of a/0 = 1 [a = /? = 3.0 x 10" 11 
cVm] and the pure Dresselhaus (Rashba) case [/3(a) = 
4.2 x 10" 11 eVm and a{(3) = 0]. The case of equal in- 
teraction strengths, when the Rashba and Dresselhaus 
interactions can cancel each other, is distinguished by 
various si gnifi cant effects reported in the literature (see, 
e.g., Ref s . TT3II25I427I) . In this case, one has the 2D electron 
system with two uncoupled spin components (see Fig. [5]) , 
each of which demonstrates the properties peculiar to a 
2DEG without the SOlJ£ Our results^ on the inelastic 
decay rate at a/ f3 — 1 arc shown in Fig. HI The sharp 
edges of the plasmon contribution are evidence of the fact 
that there is no modifications of the Landau damping re- 
gion induced by the SOI. As is seen from the inset of the 
figure, due to the shifting property E^ + = -E^Q+k- , where 
Qx = Qy = yf&m*ot, the T + and T_ curves coincide and 
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FIG. 5: (Color online) (a) Inelastic decay rate T s as a func- 
tion of k at several values of the polar angle (p^ in the small 
electron density case at a/f3 — 3.7. Dotted lines represent 
r s without the plasmon contribution (see the text), (b) The 
corresponding energy bands Ek s (at the left) and the same 
F s as a function of E^ a measured from the Fermi energy (at 
the right), (c) Spin orientations in the momentum plane for 
the considered electron gas. 



have the form of that in a 2DEG without the SOI. 

In the pure Rashba or pure Dresselhaus cases (see 
Fig. QI, the resulting r s does not tell the difference be- 
tween spin orientations in the momentum plane, which 
correspond to the Rashba or Dresselhaus SOI (sec Fig. [5]). 
As well as before, we have the relative shift (but angle- 
independent) on the momentum scale and main modifi- 
cations induced by the SOI in the energy region, where 
a quasiparticle can decay into plasmons. 

All the considered cases meet the condition of En n o <C 
Ep , where EunD = m*(\a\ + |/3|) 2 /2 is the measure of 
influence of the SOI on the band structure. However, as 
is partly discussed in Ref. [l(| in two-dimensional elec- 
tron systems with much greater En n ]j as compared to 
Ep the inelastic decay rate can substantially differ from 
that in the 2DEG without the SOI. A striking example of 
such a system is that formed by surface-state electrons in 
ordered surface alloys^ which are very promising mate- 
rials for spintronics applications. In order to predict how 
the inelastic decay rate can behave in a system, where 
Erud ~ Ep, we consider the hypothetical case with the 
unchanged m* = 0.023, e = 14.55, and /3 = 1.6 x 10" 11 
cVm, but with a = 5.9 x 10" 11 eVm (a/ (3 = 3.7) and 
n 2 D = 2.04 x 10 10 cm" 2 , which give Ep = 1.0 mcV and 
EruD = 0.85 meV. The obtained results are presented in 
Fig. [5] 

The main feature we would like to note first is that 
for holes the decay rate T_ as a function of k has an 
"outgrowth" at k < k min , where k min = m* [a 2 + (3 2 + 
2aj3 sin(2cpk)] 1//2 — the momentum, at which i?k- has a 
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minimum. A close analysis of the imaginary part of the 
screened interaction W° and the region of integration in 
Eq. (|7|) has shown that the outgrowth is caused by open- 
ing of the plasmon decay channel for transitions between 
the s = — and s' = + subbands. It is important that 
in a 2DEG without the SOI such a channel is impossible 
for holes i2i 

As is evident from the figure, in this case on the en- 
ergy scale we have a strong anisotropy of the inelastic 
decay rate. Also the presented curves clearly demon- 
strate that the latter depends strongly on the subband 
index s of the spin-orbit split band. Keeping in mind 
that the index s distinguishes spin components, we can 
say that the subband-index dependence reflects a spin 
asymmetry of the inelastic decay rate r_/r+ for a given 
direction. The asymmetry shows its worth most brightly 
in the y>k = 37r/4 direction. In fact, in that very direc- 
tion there are significant distinctions in T_ and r + as 
functions of the exciting energy and, as a consequence, 
in the corresponding IMFP for electrons. For instance, 
the ratio r_/T+ is about 3 at ~ 0.5 meV and about 2 
at ~ 1.0 meV. At further increasing of energy, the ratio 
continues to decrease. 

Note that the IMFP spin asymmetry makes a basis 
of the spin filter effect observed in hot electron transport 
through a ferromagnetic (see, e.g., Refs.0 and[13). In the 
considered case of the quantum well, the spin asymmetry 
is not such big as in ferromagnetics (see, e.g., Ref. l3l| ). 
but, as distinct from the latters, values of the IMFP spin 
symmetry depend strongly on direction and can be tuned 
by external electric field. 

IV. CONCLUSIONS 

In conclusion, we have presented a study of the in- 
elastic decay rate of quasiparticles in a two-dimensional 
electron gas with the fc-linear spin-orbit interaction that 
includes both Rashba (interaction strength a) and Dres- 
selhaus (interaction strength (3) contributions. In this 
study, the electron gas is characterized by material pa- 
rameters suitable for [001]-grown InAs quantum wells. 
We have considered the cases of a > /3 > 0, a = /3, a ~ 



(/? > 0), and (3 = (a > 0). The cases meet the condi- 
tion of E RnD < E F , where E RnD = m*(\a\ + \/3\) 2 /2 is 
the measure of influence of the spin-orbit interaction on 
the band structure. 

As compared to a two-dimensional electron gas with- 
out the spin-orbit interaction, we have revealed a relative 
shift of the inelastic decay rates for different subbands of 
the spin-orbit split band on the momentum scale. Also, 
except for the case of equal interaction strengths, we have 
found a some smoothing of sharp forms of the peak con- 
cerned with opening of the plasmon decay channel for 
electrons. We have shown that, on the energy scale, in 
this very region distinctions between the decay rates for 
different subbands become noticeable. These distinctions 
depend on the polar angle tpk and cause the inelastic 
mean free path to be angle- and subband-dependent. As 
to the case of a = f$, due to the shifting property, the 
decay rate as a function of energy has the form of that 
in a two-dimensional electron gas without the spin-orbit 
interaction. 

In order to predict how the inelastic decay rate can 
behave in a system, where E Rn o ~ Ep, we have con- 
sidered the hypothetical case of small electron density. 
We have revealed that in such a system the decay rate 
demonstrates strong anisotropy and subband dependence 
within all the considered interval of momenta and excit- 
ing energies. Since the subband dependence can be inter- 
preted as a spin asymmetry of the decay rate in a given 
direction of k, one can expect the spin-filter effect driven 
by externally applied electric field. Also, we have found 
that in the system with E Rn p> ~ Ep holes can decay 
into plasmons, what is impossible in a two-dimensional 
electron gas without the spin-orbit interaction. 
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